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We calculate the Kondo temperature and crystal-field levels of strongly correlated multior- 

bital systems solving the Anderson Impurity Model with the finite U Non-Crossing Approximation 
(UNCA) in its simplest scheme, that is, considering the self energies at lowest order in the 1/N 
diagrammatic expansion. We introduced an approximation to the vertex function that includes the 
double energy dependence and investigate its effect on the values of Tk for simple electronic models. 
We also analyze the competition between the two spin fiip mechanisms, involving virtual transitions 
to empty and doubly occupied states, in the determination of the ground state symmetry by includ- 
ing an extra diagram of higher order in 1/A''. We finally combine the resulting simple formalism with 
ah initio calculated electronic structures to obtain Tk's, ground states, and crystal field splittings 
in excellent agreement with experimental results for two particular Ce compounds, namely Celna 
and CeSna. 



I. INTRODUCTION 

One of the most studied models for strongly correlated 
electron systems is the Anderson impurity model. This 
can account for the Kondo regime, where conduction elec- 
trons scatter off a localized magnetic impurity and form 
a local singlet. Initially proposed as a model for local 
moment formation in metals, in recent years it has been 
extensively used to describe the physics of quantum dots 
and strongly correlated lattice systems, mainly, heavy 
fermion compounds. 

Among the several methods proposed to solve the An- 
derson impurity Hamiltonian, the so-called non-crossing 
approximation (NCA), in its lowest order self-consistent 
form occupies a special place due to its computational 
simpHcity. NCA can be thought as a pertubative ex- 
pansion with respect to where TV is the degener- 
acy of the impurity levels. In particular, the NCA has 
been widely used in the last twenty years to solve the 
Anderson Hamiltonian in the infinite U limit, in which 
the double occupancy of the impurity site is prohibited. 
This scheme was successfully applied to, both, theoreti- 
cal models and real materials jl,|5i,i6i], to analyse 
different physical properties, like magnetic susceptibili- 
ties, crystal-field splittings and spectral properties. Also, 
within this approximation, the out-of-equilibrium trans- 
port properties of single and double quantum dot systems 
have been studied[7, 8]. 

For U ^ CO and large N degeneracy of the local states, 
NCA captures the Kondo energy scale and provides a 
qualitative description of the formation of the Kondo 
resonance when the temperature approaches the Kondo 
regime [T Tk) from above However, this is not 
the case for temperatures much lower than Tk {T <C Tk), 
where NCA yields unphysical results. Furthermore, the 
Tr-'s are often underestimated in the infinite U limit. 

The deficiencies that restrict the usefulness of NCA 



are greatly reduced if a finite on-site repulsion U is con- 
sidered together with the inclusion of vertex corrections 
Q. While NCA in the infinite U limit contains all non- 
crossing diagrams up to the (l/A^)^^^ order, with a finite 
U there appear crossing diagrams of order (l/iV)(o) which 
have to be at least included through vertex corrections. 

Recently, Sakai et al. [l3| and J. Otsuki and Y. Ku- 
ramoto have shown that a finite-C/ NCA (UNCA) 
that collects all the order diagrams yields the 

proper energy scale of the Kondo effect. Within their ap- 
proximation, which also considers non-crossing diagrams 
of order the main problem is how to achieve a 

complete solution of the integral equations for the ver- 
tex corrections, that depend on two independent energy 
variables. 

A more sofisticated extension of the NCA technique 
has been introduced by Haule et aZ.[l3l- This exten- 
sion, known as the symmetrized finite-?/ NCA (SUNCA), 
treats the fiuctuation processes into the empty and into 
the doubly occupied intermediate states on an equal foot- 
ing, by means of a proper symmetrization of the vertex 
corrections. Although SUNCA provides a correct energy 
scale, its practical computation is not easy, even for a 
simple structure of the conduction electron band. 

In order to be applied to real materials, the scheme to 
solve the many-body Hamiltonian should be able to give 
both, a good approximation to the Kondo temperature as 
well as the correct ground state symmetry without imply- 
ing an unreasonable computational effort under realistic 
conditions (multiorbital models, complex band electron 
structures, etc.). 

In an earlier contribution, we have been able to under- 
stand and predict trends in the evolution of crystal-field 
splittings by using the NCA in the infinite U limit, with 
the hybridization function taken from ah initio electronic 
structure calculations The aim of this contribution is 
twofold, on the one side we introduce an approximation 
to vertex functions including the double energy depen- 
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dence in two different ways within UNCA at lowest or- 
der in and study its influence on Tk- On the other 
hand, we analyze the importance of diagrams of differ- 
ent orders on the correct prediction of the ground state 
symmetry. Finally, we exemplify with two real Cerium 
compounds whose hybridization functions are obtained 
from flrst principles. 

The paper is organized as follows. In section II, we 
introduce the Anderson Hamiltonian and the auxiliary 
particle formahsm. We also present the UNCA equations 
up to the order of the diagrams. In section III 

we summarize the basic results for the Kondo temper- 
ature implementing different treatments for the vertex 
functions. In section IV, we investigate the crystal-fleld 
splittings induced by hybridization and the symmetry of 
the ground state. Finally, in section V, we apply the pre- 
vious ideas on two real systems. The Kondo temperature 
and the symmetry of the ground state are obtained for 
Celna and CeSns. In section VI, we summarize and con- 
clude. Some details of the calculations can be found in 
the Appendices. 



|m) = 4» \vac) , 

where \vac > is the vacuum state for the auxiliary parti- 
cle operators, b'' is the "light" boson (empty state), sj„'s 
are the single pseudo-fermions (single-occupied states 
with energies Em), and rfm,„/'s are the "heavy" bosons 
(which correspond to the doubly occupied states with 
energies E„in — f-m + £n + U) . These heavy bosons sat- 
isfy the antisymmetric property d^^^, ~ ^d^m'm- '^^^ 
local physical electron operator, /„, can be written as 

,/m b^ ^ni ^" ^ ^ ^ffi'^jnm' ■ 

In this representation the Hamiltonian becomes 

km m 



II. AUXILIARY PARTICLE REPRESENTATION 
AND THE LARGE-Af NCA 

The Anderson impurity model with finite U is de- 
scribed by the Hamiltonian 

H — efcm (^km^km 

km 

m m<n 

+ Y [VkrafLckm + Vfc*„cj,„/„) , (1) 
km 

where the indices to, n label the quantum numbers of 
the impurity levels. The operators c^^, create a con- 
duction and a localized electron state, respectively, and 
^m — fmfm is the /-number operator. The last term 
represents the hybridization between conduction and lo- 
cahzed electrons, and Vkm are the hybridization matrix 
elements. In the auxiliary particle approach to the An- 
derson impurity model, the local impurity states are 
represented by additional degrees of freedom, assumed 
to be created by pseudo-bosons and pseudo-fermions 
operators Q. In order to take into account the effects 
of a finite value of the on-site Coulomb interaction, U, 
the local Hilbert subspace must contain the empty, sin- 
gle, and doubly occupied states, while states with higher 
occupancies can surely be neglected for intermediate to 
large values of U. The local states are represented as fol- 
lows: 

|0) = fot \yac) , 



+ ^ (Cm + + U) dl^„dmn 

+ Y {^kmslnb Ckm + H.C.^ 
km 

+ ^ ^ (^kmdlnm'Sm'Ckm 

+ H.C.). (2) 

kmm' (m' ^m) 

The auxiliary particle Hamiltonian is invariant under 
simultaneous, local U{1) gauge transformations, a — > 
de"^^*-*, where (pit) is an arbitrary, time-dependent phase 
and d is any auxiliary particle operator. Due to this, 
there is a conservation of local charge Q in time. There- 
fore, all physical quantities should be obtained with the 
constraint Q = 1 0]. The corresponding Green's func- 
tions of the auxiliary particles have the usual structure 
Ga^i^) = z — Ca — T,a{z), and the self-energies T,a{z) can 
be evaluated by means of the non-crossing approxima- 
tion. 

In the context of the NCA, the simplest picture 
that contains the two elementary "spin-fiip" scattering 
processes, involving empty and doubly occupied in- 
termediate states, respectively, and that captures the 
correct energy scale of the Kondo effect in the finite U 
case, is obtained with the inclusion of vertex corrections 
in a large- TV expansion, retaining only the lowest order 
diagrams in (1/A^), that is, (l/iV)(°).[l3|. In this approx- 
imation, called from now on UNCA(°), the self-energy 
of the heavy boson propagator vanishes. Consequently, 
the Green's function of the doubly occupied states can 
be written in the simple form Gmln'i^) = z — E^m' ■ On 
the other hand, the self-energy for the pseudo-fermion 
propagators contain only the contribution coming from 
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Figure 1: Throughout this article, the full, wiggly, dashed, 
and curly lines stand for pseudofermion, empty boson, con- 
duction electron, and heavy boson propagators, respectively. 



+ 



= I + 



Figure 2: Dyson's equation for the pseudo-fermion and light 
boson propagators within the non-crossing approximation up 
to (1/A'')'' order. The lines used in Feynman diagrams are 
shown in Fig. [T]The bold (light) lines represent the full (free) 
propagators. The big dot represents the vertex function, that 
is, the bare hybridization vertex Vkm dressed by vertex cor- 
rections. 



equations: 



— nFie)Tm'{e)Gmm'i^ + e), (3) 



/de' 



X Gm' -t- e')Smm' + (-' + 



(4) 



/d€. 
— npie) T„i{€)Am{i^,e)g.m{uj + e). (5) 



Here, npie) is the Fermi function and Tm{^) are the 
hybridization functions between the conduction electron 
band and the impurity state. 



(6) 



It is easy to see that, with the self-energies for the pseudo- 
fermions computed from Eq. ([3]), one can obtain the 
vertex functions through Eq. ^ . Finally, the self-energy 
for the light boson propagator is obtained directly from 
Eq. A complete solution of the equations for the 

vertex functions seems to be impractical. In the next 
section we introduce some simplifications in order to solve 
them numerically with low computational effort, even for 
realistic systems. 



III. APPROXIMATIONS TO THE VERTEX 
FUNCTION 




Figure 3: Diagrammatic representation of the integral equa- 
tion for the vertex function in the ladder approximation. 
The vertex corrections incorporate crossing diagrams of or- 
der {l/Nf to the NCA equations. 



the heavy boson. The crucial advantage of this scheme is 
that the system of equations does not have to be solved 
in a self-consistent way. 

The self-energies, denoted by loops in the diagrams of 
Fig. [21 can be evaluated by standard Feynman rules and 
by the projection process into the Q = 1 subspacefl3|- 
The vertex function Am(ti;,e), in the ladder approxima- 
tion shown in Fig. O can also be evaluated in the same 
way. The final result is the following set of non-coupled 



The vertex functions depend on two independent en- 
ergy variables, uj and e. A complete solution of the set 
of integral equations @ implies a lot of computational 
effort, hence some simplification is needed. The schemes 
NCA and UNCA(°) are conserving, that is, the diagrams 
can be derived from a Kadanoff-Baym functionaljla. flB|. 
However, this property is lost when approximating the 
vertex function. This non-conserving feature, does not 
imply major consequences if the spectral functions of the 
slave particles are always positive and also integrate to 
one, as expected for the conservation of the expectation 
value of the number of auxiliary particles in this theory 
(Q). We check that this is the case in the following cal- 
culations. 

In this section we compare the simplest approx- 
imations to solve numerically the vertex equations 
and analyse their infiuence on the calculated Kondo 
temperature for different values of the on-site Coulomb 
repulsion U . 
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A. NCA/^v approximation 

Sakai et aZ[l3| have pointed out that the exchange cou- 
phng due to virtual transitions to the doubly occupied 
state can be obtained even when the e energy depen- 
dence of the vertex function A™ (a;, e) is neglected. So, in 
order to get rid of the variable e, they replace the heavy 
boson propagator in Eq. Q by a representative value, 
that is, Gmm'{z ~ Em') = — l/(em + U) ■ Within this ap- 
proximation there is no more e dependence in the vertex 
functions: 



In this way, Eq. (U is simplified as follows 



where 



de 



U 



Fm{uo) = — nF(e)r™(e)C?„i([j -I- e). 

' TT 



(7) 



(8) 



(9) 



With this simplification of the vertex functions, Am\ the 
boson self-energy turns into the simple expression 



(10) 



equations ([4]) become a set of linear algebraic equations 
for given values of u and e, 

J (12) 
Tm' (e')5m' (w -I- e')^mm' {uj + e' + e) . 

To further simplify the integral, we can evaluate the dou- 
bly occupied Green's functions at some specific value of 
Lu, as before. The energy dependence on e for the dou- 
bly occupied Green's functions could be kept if we set 

LO + e' = e„i/ , 



Gmm' (uj + e' + e) 



-1 



U — e — irj 



The final form of the vertex function is again obtained 
from a set of linear algebraic equations. This and the cor- 
responding self-energy of the boson propagator are given, 
respectively, by 



and 



(13) 



(14) 



Sakai et al. have called their approximation NCA/^v 
and, using it, they have obtained good results for 
different properties of real systems, like CeSb^Ol]. On 
the other hand, Svane et al.\v^ have also used vertex 
functions with a single energy argument to analyse the 
monopnictides series of Cerium compounds, CeX (X=N, 
P, As, Sb, Bi). 



To be consistent with the nomenclature introduced 
by Sakai, we call this approximation the 'NCA/^ve 
approximation', where e indicates the double energy 
dependence of the vertex functions. 



C. NCA/ ve approximation 



B. NCA/^ve approximation 

We can recover the second energy dependence, e, in 
the vertex functions, avoiding at the same time to solve 
the integral equations, if we interchange energy variables 
in the integrand of Eq. ^ as follows: 

A^,{uj,e')^Al^]{co,e). (11) 

This procedure is justified under the assumption that the 
e dependence of the vertex functions is relatively weak. 
The resulting approximation for the vertex functions was 
successfully used by Kang et aZ.[JJ;] to study the infiuence 
of a magnetic field in the impurity Anderson model. Even 
though we recover the double energy dependence of the 
vertex functions, it should be noticed that energy conser- 
vation is still not completely satisfied in the vertex cor- 
rection diagrams. Under this approximation, the vertex 



The NCA/^ve vertex function can be further improved 
by evaluating the doubly occupied Green's functions (Eq. 
(fT2| ) at the value of w -I- e' that most likely maximize 
the rest of the integrand, that is, uj + e' ^ em', where 
im' are the renormalized poles of the pseudo-fermion 
Green's functions. From Eq. ([3]), we compute directly 
the imaginary parts of the pseudo-fermion propagators, 
that are the spectral functions corresponding to the re- 
tarded Green's functions Gmiz) 

Pm(^) = --Im g„i{oj). 

TT 

The poles of Pm{uj) are the energies im that satisfy 

£m ~ ^ Sm(e„i) = 0. (15) 

If we set ijj + e! — in the argument of the doubly occu- 
pied Green's functions in Eq. (fT2l) . the vertex functions 
and the boson self-energy are then given, respectively, by 
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A(^)(a;,e) + ^— i 



(16) 



(17) 



D. Kondo temperature results 

We analyze next the effect of the different approxima- 
tions described above on the Kondo temperature, within 
a simple model. For this purpose, we use a constant 
and degenerate hybrization intensity Tm{() — 0.15 eV 
for —B < e < B and otherwise. Here B is the half 
bandwidth and we set S = 3 eV. The degeneracy N in 
this section is taken to be = 6 and we set e™ — —2 eV 
for all m's. With these parameters we solve the UNCA(*" 
set of equations in the NCA/^v, NCA/^ve, and NCA/^ve 
vertex approximations. We obtain the spectral functions 
for the boson (p6(w, T)) and pseudo-fermions (pm(i^, T)), 
in the T ^ limit, in order to calculate the Kondo tem- 
perature, Tk- As usual, this temperature is obtained 
from the difference between the lowest pole of Pb(w) and 
the corresponding one for Pm('^)[l!|- We obtain Tk for 
different values of the Coulomb interaction U and within 
the different approaches to the vertex corrections previ- 
ously introduced, taking U= 5 eV, 10 eV and 100 eV. 
We consider that U= 100 eV already gives the U oo 
limit. The U ^ oo limit is characterized by Am — 1 (to 
leading order of the large- A^ expansion) and this limit 
is recovered perfectly with U = 100 eV. Our results are 
presented in Table [H 



U 



r, 



(a) 



(6) 



5 

10 
100 



260 
85 
33 



208 
81 
34 



183 
77 
31 



Table I: Results for the Kondo temperature obtained from 
different approximations to the vertex functions. The super- 
scripts a, b and c represents the NCA/^v, NCA/^ve, and 
NCA/^ve approximations to the vertex functions respectively. 
The Kondo temperatures are given in Kelvin. The values of 
U are given in eV. 

In Fig. (jU we compare the real parts of the vertex 
functions in the NCA/^v and NCA/^ve approximations 




w(eV) 

Figure 4: Real part of the vertex functions as a function of uJ 
for f/ = 5 eV. The second argument of Am\uj, e) is evaluated 
at the pole of the boson spectral function, Eq. 



as a function of to. To plot the function i?eAm''(i^, e), we 
set its second argument equal to the value of the lowest 
pole of the boson spectral function, Eq. The correspond- 
ing curve for the vertex function in the NCA/^ve ap- 
proximation, aI^{lj,Eq), is very similar to Am\uj,Eo), 
so that it is not shown in Fig. From this figure it 

can be seen that the NCA/^v vertex agrees qualitatively 
with the NCA/^ve one. As it can be drawn from TablelH 
all the approximations give the same order of magnitude 
for the Kondo scale, for a given U. However, there is a 
quantitative difference in the calculated Tk 's that can go 
up to 30 %. 

On the other hand, in Fig. ([5]) we plot Am{Eo, e) and 

A^{Eo,e) as functions of their second arguments. For 
u! we have chosen, just as an example the value of the 
boson propagator pole, Eq. Both curves seem to have 
been rigidly shifted one with respect to the other, as a 
consequence of the different argument values in which 
the Green's functions of the heavy bosons are evaluated 
in the vertex corrections. The differences show up in the 
values obtained for Tk- 

Comparing Figs. (0]) and |[5]), it can be observed that, 
as a function of w, there is an important variation of the 
functions in a wide energy range while, as a function of 
e, the energy window of variation is narrower. However, 
in this last case the amplitude of variation is up to five 
times greater than in the first one. From this fact, it 
is clear that the second dependence in energy cannot be 
neglected for J7 = 5 eV if precise values of Tk are desired. 
This result is relevant in the case of Ce systems for which 
[/ - 6 eV. 
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Figure 5: Real part of the vertex functions in the NCA/^ve 
and NCA/^ve approximations as a function of e for (7 = 5 
eV. The first argument of both functions is evaluated at the 
pole of the boson spectral function, Eq. 



IV. CRYSTAL-FIELD EFFECTS 
A. Model 

In many real systems the crystal-field splittings are 
important and have not been taken into account in the 
previous model. In this section, we investigate a model in 
which the hybridization functions r„i(e), defined by Eq. 
([6]), are energy independent, but have different intensities 
according to the corresponding m-symmetry. Therefore, 
in these cases there are hybridization-induced crystalline 
electric field splittings, A„„/. The splitting induced by 
hybridization is anticipated to be the dominant contribu- 
tion to the crystal-field splittings in strongly hybridized 
Kondo systems, such as Cerium compounds|l8|. The 
crystal-field splittings are defined as the absolute value 
of the difference between the lowest pole of the pseudo- 
fermion spectral functions (e„) and the excited ones (e„'), 
A„„' = |en' — Enl- We consider the A„„/ coming from vir- 
tual transitions /" ^ and ^ ^ /i of the 
local level, that is, we focus on the splitting coming from 
the anisotropy (m-dependence) of the mixing interaction 

A good scheme to solve the Hamiltonian should be able 
to give, not only the correct Kondo energy scale but also 
the correct symmetry of the local impurity ground state. 

We can get a preliminary idea about the symmetry of 
the ground state as shown in Ref. Within the infi- 

nite U NCA[^, the pseudo-fermions are dressed only by 
the empty boson. For a fiat density of states of the con- 
duction electrons, of half- width -B, the corresponding real 
part of the self-energy of the pseudo-fermions is given, at 
zero temperature and to leading non-trivial order 1/iV, 



by the approximate relation 



U=oa,NCA 



' 7 

-In 



(18) 



with K a positive real number. Therefore, the shift of the 
pseudo-fermion energies is always negative and directly 
proportional to the hybridization strength. Using this 
simple argument, we can foresee that the level with the 
strongest hybridization has the largest shift and becomes 
the ground state. It is worth to notice that the same 
physics holds in the resonant model (C/ = 0). 

On the other hand, within UNCA^") , the shift of the 
different pseudo-fermion levels is given by the self-energy 
coming solely from the doubly occupied bosons (see Eq. 
([3])). In order to compare the A„„''s coming from empty 
or doubly occupied bosons, we estimate this last con- 
tribution. We analyse here a system with two different 
constant hybridizations, namely Fi and which are A^i 
and N2 times degenerate, respectively. The real part of 
the self-energies of the pseudo-fermions are given by the 
following approximate relation, at zero temperature. 



U<oo 



E ■ 

m' 



u 



U + B 



(19) 

where a is also a positive real number and considering 
that all the have the same value. We have found, 
taking into account only virtual transitions — > /"^ ^ 
{U ->oo limit). 



U=oo,NCA 



ReY, 



U=oo,NCA 



-kTi 



-KF2, 



(20) 



but on the other hand, taking into account only virtual 
transitions f^^p-^f^,{U< 00), the self energies are 
given by 



U<oo 



-ai{Ni-l)T,+N2T2) 



-aiNiTi + iN2-l)T2 



(21) 



Considering the case Fi > F2, the corresponding self- 
energies are related to each other as follows 



n ^U=oc,NCA ^ r> ^U=oo,NCA 



U<oo 



> ReY, 



U<oo 



(22) 



(23) 



Therefore, we can infer that there exists a competition 
between both spin fiip mechanisms, ^ Z*' — > and 
P f^, to give the ground state symmetry, as it 
will be confirmed in the next subsection. 

Due to this, we add to the UNCA(°) dia grams the 
fermion self-energy diagram of order (1/iV)^, originated 
in virtual transitions to the empty state, and we analyze 



7 



+ 



+ 



+ 



Figure 6: Diagrammatic representation of the finite U non- 
crossing approximation up to {I/NY' order with an additional 
diagram of order {l/Ny. The big dot represents the vertex 
function. The diagram of order (1/A'')^ has a bare light boson 
propagator, and it does not contain vertex corrections. 



the influence of this last contribution on the A„„''s values 
and on the ground state symmetry. 

In Fig. ^ we show the diagrammatic representation of 
the non-crossing approximation up to (UNCA'-^-') 
order with the additional diagram of order (1/iV)^. We 
call this approach UNCA(*'+^^ from now on. In order to 
compare on an equal footing both contributions to the 
self-energy of the pseudo fermions, in this additional dia- 
gram we neither consider vertex corrections nor the self- 
consistent renormalization of the light boson propagator. 
Summarizing, in this approximation, the self-energy of 
the pseudo-fermions has two contributions, one coming 
from the empty boson and another from the doubly oc- 
cupied ones, as shown by the following expression 



E 



de 



n_F(e)rm'(e)^mm' -l- e) 



— [l-nF(e)] 

TT uj — e 



B. Crystal-field results 



(24) 



As a flrst step in the modeling of real compounds, 
in this subsection we consider a system with two differ- 
ent constant hybridizations, Fi and F2, with degeneracy 
Ni = 2 and A'2 = 4, respectively. As before, we take 
the half bandwidth B = 3 eV and £,„ = —2 eV for all 
m's. To obtain the Ai2's we compute the pseudo-fermion 
self-energy coming from Eqs. ^ and l(24|) . cor- 

responding to UNCA^") and UNCA(°+i), respectively, 
which are solved for different values of U. In this subsec- 
tion we use the vertex corrections given by the NCA/^ue 
approximation. 

In Table |ll] we show the calculated values of A12 for 



U 



A !7iVCA(°' 
^12 



^12 



5 

10 
100 



24 
12 
2 



6 
21 
31 



Table II: Results for the crystal-field splitting A12 = |ei — £2! 
obtained for AF = 0.1 eV within the NCA/^ve approximation 
to the vertex functions. The values of the A12 are given in 
Kelvin. The values of U are given in eV. 



Fi = 0.15 eV and F2 = 0.05 eV. We obtain the F2 sym- 
metry as the ground state when we use UNCA^"-* and a 
decrease of the absolute value of A12 when going from 
[/ = 5 eV to [/ ^ 00. It is clear from Eq. that 
a vanishes in the U ^ 00 hmit and that the shift pro- 
duced by S^^°° becomes equal to the one corresponding 
to S]^"^°°. The crystal-fleld splittings in the inflnite U 
limit come only from the ^ f'^ ^ processes, as 
expected. 

On the other hand, UNCA(°+^) gives Fi as the ground 
state symmetry, as it has been anticipated by the ap- 
proximate expression of Eq. (|23| . In this case, we obtain 
an increase in the values of A12 when going from U = 5 
to U -^00. Similar behaviour is obtained for different 
values of AF = Fi — F2. 




10 15 20 25 30 35 40 45 50 
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Figure 7: Kondo temperature, Tk, as a function of the crystal- 
field splitting, A12, within the UNCA'°' scheme. The ob- 
tained symmetry for the ground state is r2 in all calculations. 
The values are given in Kelvin. 



We then analyse the evolution of for different values 
of AF by setting Fi = 0.15 eV and varying F2, from 0.09 
to 0.15 eV. 
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Figure 8: Kondo temperature, Tk, as a funcions of the crystal- 
field splittings, Ai2, within the UNCA'°+^' scheme. The ob- 
tained symmetry for the ground state is Fi in all calculations. 
The values are given in Kelvin. 



For each AF we calculate the splitting A12 and Tk 
within UNCA(°) and UNCA(°+i) as shown in Figs. © 
and ([8]) respectively. In both cases, and for all the val- 
ues of U considered, we obtain a decrease of the Kondo 
temperature as the crystal field splitting increases. In 
this simple context using constant hybridization func- 
tions, the increment in the splitting of the levels reduces 
the effective degeneracy of the ground state which, in 
turn, is correlated with a decreasing Tk- This behaviour 
is more pronounced when using 

On the other hand, it has to be noticed that in the 
UNCA(°) calculation, when U is large, a small variation 
in AF gives rise to large variations in A12. This leads to 
low Kondo temperatures and even to the disappearence 
of the Kondo regime {Tk < 0) which can be clearly seen 
in Fig. lO for U= 100 eV, where only two points present 
a positive Tk- 

It must be stressed that the evolution of Tk with A12 
shown in Fig. Q and ([HI is not necessarily the one of a 
real compound. We will show in the next section that 
using ab initio calculated hybridization functions, which 
show a rich structure as a function of energy, can give 
rise to a different behavior. 



V. REALISTIC CALCULATIONS FOR CERIUM 
COMPOUNDS 

In this section we apply the previous ideas to two re- 
alistic Cerium compounds with different electronic prop- 
erties. In the case of Ce intermetaUic compounds long 
experience in the area shows that the behavior of differ- 
ent properties, such as Tk, crystal field splittings, mag- 



netic susceptibilities, can be readily obtained, qualita- 
tively as well as quantitatively, by considering an impu- 
rity problem due to the locaHzation of the 4f orbitals 
[1, @, [3, [l^, [2l|. In particular, we choose the heavy 
fermion system Celns and the intermediate valence one 
CeSna. In both systems, from mag netic susceptibility 
measurements, Pedrazzini et al [2a| have obtained the 
J = ^ Tf doublet as the ground state. The same result 
was obtained from neutron scattering measurements by 
Murani et al [131 ■ Celna shows the normal Curie- Weiss 
susceptibility at all temperatures except around and be- 
low the Neel temperature, Ti\[{= 10.2K), so that a small 
value for the Kondo temperature is expected {Tk < TV) 
[i^ l- CeSns, on the other hand, shows an enhanced Pauli 
susceptibility at low temperatures followed by a shal- 
low maximum and a Curie- Weiss behaviour above it [23 |. 
The Kondo temperature in this case is well reported and 
its value is Tk = 450X[i|. 

In both compounds a precise experimental estimation 
of the crystal-field splitting is difficult. In the case of 
Celna, In is a strong neutron absorber and this makes the 
interpretations of the experiments rather cumbersome. 
In spite of this, the reported estimated value for Acf is 
around 130 K (2^ . [23 |. On the other hand, in the case 
of CeSns the strong hybridization induces a large Tk so 
that Tk » Acf and the experiments fail to resolve the 
crystal-field peaks. 



A. Ab initio hybridization functions 

We obtain the hybridization functions of the 4f states 
of Cerium with the conduction band from first principles 
within the Density Functional Theory. In this work, the 
ab initio calculations are done using the full potential 
linearized augmented plane waves method (FP-LAPW), 
as implemented in the Wien2k code [131 ■ As suggested 
by Gunnarsson et al [IBl, the hybridization Fm(e) can 
be estimated from the projected LDA 4f density matrix 
P 



LDA 



at the Ce site in the following way, 



Tmle) = -Im < lim 



^LDA 



dz- 



(^) 



(25) 



In all cases the labels m correspond to the different 
irreducible representations of the 4f states at the cubic 
Ce site. That is, for J — ^ the doublet r7 and the quar- 
tet Fg, while for J — ^ the doublets Fg and F7, and 
the quartet Ts- The LDA calculations are performed at 
the experimental volumes of the Celns and CeSns com- 
pounds. The muffin-tin radii, Rmt are taken equal to 
2.4 a.u. in the case of the anion Hgands, while the corre- 
sponding radii for Ce are taken equal to 3.0 a.u. in the 
Celna compound and 3.3 a.u. in the CeSna one. 102 k 
points in the irreducible Brillouin zone are considered to 
be enough for the quantities to be calculated. 

The hybridization function is used as input in the 
UNCA set of equations. The crystal-field splittings are 
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read from the separation of the peaks of the different 
spectral functions, Pm's, which are shifted one with re- 
spect to the other due to the different degree of hybridiza- 
tion of each 4/ level with the conduction band. We focus 
on the value of the splitting in the = | multiplet, 
namely Ac_f = £/r7 — e/rs- We will call this mixed tech- 
nique LDA-UNCA from now on. In the UNCA equa- 
tions we take the bare energy value for the 4/ state from 
photoemission experiments [2(3|. It is, namely, -2 eV for 
Celna and CeSns. We shift the = | multiplet in an 
amount given by Aso = 0.35 eV, due to the spin-orbit 
interaction, and use U =6 eV for the on-site Coulomb 
interaction constant among the /-electrons, which is con- 
sidered to be a standard value for Cerium systems [27| . 
All these energy levels are given with respect to the Fermi 
energy. 
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Figure 9: Hybridizations functions from LDA calculations for 
the J = 5/2 multiplet of Celna. 



B. LDA-UNCA results 

The calculated hybridization functions for the J — 5/2 
multiplet are shown in Figs. ^ and ifTOj) in the cases 
of Celna and CeSns, respectively. These hybridization 
functions are in qualitative agreement with the ones cal- 
culated by Han et al. U| using the LMTO-ASA approx- 
imation. In both compounds, the r7 symmetry has the 
largest average value in the whole range of energy. Using 
the LDA-UNCA method, we obtain a ground state of Fg 
symmetry for both systems when we use the UNCAt") 
scheme. This does not agree with the experimental re- 
sults for these Cerium compounds. The Ac_f's do not 
depend on the approximation employed to compute the 
vertex functions and their values are Acf = 1%K and 
Ac_F — SOX for Celna and CeSna, respectively. On 
the other hand, within the UNCA'^°+^) scheme, we ob- 
tain a ground state with Fy symmetry for both studied 




£(eV) 



Figure 10: Hybridizations functions from LDA calculations 
for the J = 5/2 multiplet of CeSna. 



Cerium systems, and for all vertex corrections used. In 
this scheme, the values of the Acf are 150 K and 180 K 
for Celns and CeSna, respectively. In the first case this is 
in good agreement with the experimental result reported 
for Celna (130 K). As we have shown in Section III D, 
the Tk strongly depend on the choice of the approxima- 
tion for the vertex functions. For these real systems, the 
Tk are again overestimated in the NCA/^v approxima- 
tion with respect to the NCA/^ve one. In Table Hill we 
show the calculated values for the Tk , as well as the sym- 
metry of the ground state, for the CeSns compound in 
both, UNCA(o) and UNCA^o+i) schemes, and for ah the 
approximated vertices considered. 

We want to remark that the UNCA^^"'"^-' scheme with 
the NCA/^ve approximation for the vertex corrections 
gives values for the ground state symmetry and for the 
Acf's in very good agreement with the experimental 
ones in these Cerium systems. Furthermore, for the 
CeSna compound, we obtain a value for Tk very near to 
the experimental one, that is,450i^r. In the case of Celna 
we obtain a negative value for the Tk which implies that 
the Kondo regime has not been completely established in 
agreement again with the experimental observation that 
Tft:, if it exists, should be below Tjv ~ 10 K as mentioned 
above. 

All the vertex approximations presented in this work 
capture the experimental energy scale of Tk in the case 
of CeSns, but only the NCA/^ve approximation within 
the UNCA(o+i) scheme gives a quantitative agreement 
{Tk = 480X) with the experimental value {Tk = AbQK). 

As already mentioned, the results obtained for Tk and 
for the crystal-field splittings in these compounds do not 
follow the behavior of the constant hybridization simple 
model of the previous section, in the sense that a larger 
splitting does not necessarily lead to a smaher Tk- This 
enforces the importance of taking into account the real 
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CeSns Tgs NCA/^v NCAfve NCAfve 
Tk Tk Tk 



UNCA'") Ts 1030 788 722 

UNCA'°+^' Ft 757 538 480 



Table III: Results for the Kondo temperature and symmetry 
of the ground state (res) as obtained using different approx- 
imations to the vertex function. Within the UNCA''^-' and 
UNCA'""''^-' schemes for the CeSns compound. The Kondo 
temperatures are given in Kelvin. 



electronic structure information when solving many-body 
hamiltonians. 



VI. DISCUSSION AND CONCLUSIONS 

We have solved the Anderson impurity Hamiltonian 
using the finite U non-crossing approximation (UNCA) 
in its simplest scheme, that is, considering the self ener- 
gies at lowest order in the 1/N diagrammatic expansion 
(UNCA(°)). This approximation yields a correct Kondo 
energy scale, and it is then a reliable approach to perform 
realistic calculations for 4f systems. 

Even if the approximation is simple, a set of coupled in- 
tegral equations for the vertex functions has to be solved. 
These functions depend on two independent energy vari- 
ables. In this work, we have extended the NCA/^v ap- 
proximation to the vertex functions introduced by Sakai 
et al, including the second energy variable in two differ- 
ent ways called NCA/^ve and NCA/^ve. We have cal- 
culated and compared the Tk using these three different 
approximations. Our results confirm that the NCA/^v 
approximation gives the correct Kondo energy scale, but, 
in order to make a quantitative comparison with the ex- 
periments (in real systems), the double dependence on 
energy cannot be disregarded when the Coulomb inter- 
action involved is near to U = 6 eV. In fact, with our 
NCA/^ve approximation, we obtain Tk for the CeSna 
compound in surprising agreement with the experimental 
one, when using ab initio calculated hybridization func- 
tions. 

We have also analyzed the competition between the 
two spin fiip mechanisms, involving virtual transitions to 
empty and doubly occupied states, in the determination 
of the ground state symmetry of the local state, by in- 
cluding an extra diagram of higher order in This 
extra diagram corresponds physically to virtual transi- 
tions to the empty state, — > /^^ — > which is not 
included in the UNCA^^^ scheme. We have shown that 
the inclusion of this diagram (UNCA^^"*"^^) is necessary 
in order to obtain the experimental ground state symme- 



try for the Cerium systems Celns and CeSns, which we 
have taken as test examples. 

Summarizing, this work presents an analysis of the dif- 
ferent corrections that can be made on the vertex func- 
tions within the finite-U NCA approximation showing, 
quantitatively, the importance of its double energy de- 
pendence on the Kondo temperature. It also gives phys- 
ical insight into the factors determining the ground state 
symmetries as a competition between two different spin 
fiip mechanisms. We finally test and apply the previ- 
ous features in real systems taking as input the LDA 
hybridizations, obtaining very good agreement with ex- 
perimental information. 

VII. APPENDIX 

A. Numerical Details 

All vertex corrections and self-energies have been com- 
puted directly from their expressions given in this paper. 
To calculate them, we employ two different sets of fre- 
quencies, e(i) and with 1 < i < A'h and 1 < j < Ah 
respectively. The mesh for e(i), in which the hybridiza- 
tion function is defined, is divided into Mi equal steps of 
size h = {cmax — emin)/Mi. We use a large value for Mi, 
typically 4000. In view of this, a simple trapezoid rule 
is enough to perform the integrals. On the other hand, 
the mesh is a non-hnear one and collects the major- 
ity of its points near the pole of the light boson spectral 
function. This mesh is also dense enough near the pole 
of the pseudo-fermions spectral functions. This attribute 
becomes fundamental to give a well defined Tk work- 
ing with a reasonable amount of w frequencies. In this 
analysis, we use a number M2 of frequencies around 
1000-1500. The infinitesimal 77 in the Green's functions 
is fixed as a multiple of the minor step in the uj{j) mesh. 

B. Linear algebraic equations for PSm{Lo,e) 

In this subsection, we give some details about the set 
of linear algebraic equations given by Eq. iflG]) . We use 
the short notation 

em(e) = em + U - e-i-q 

for the equation that defines the vertex Am\ which be- 
comes 

in' ^tn 

The functions Fm(aj) are given by Eq. ([9]). 

In the specific case in which the index m runs over 
three different symmetries, namely [1, 2, 3], whose degen- 
eracies are {2, 4, 8} , respectively, the solution of the pre- 
vious system requires the inversion of a complex non- 
symmetrical 3x3 matrix. This inversion is done for each 
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and e{i) energies. 

AW(w,e) = M-^(a;,e) e{e) 

where the matrix M, the vectors A('^^(a;,e) and £(e) are 
given by the following expressions 



M{uj,t 



This is the tj'pical case of Cerium systems in which 
level TO = 1 and to = 2 correspond to the Fy, Tg 
symmetries of j = 5/2 multiplet, respectively, and 
to = 3 represents the whole j = 7/2 multiplet. 



^ Fi{oj) + ei{e) 4^2 (w) 8^3(0;) ^ 

2Fi(a>) 3F2(a)) +c2(e) 8i^3(w) 
\ 2Fi{u) 4^2 M 7i^3H+e3(e)/ 



A(^)(w,e) 



2 ('^^e) 



A, 
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